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Abstract 

We study the XXZ chain with a boundary at massless regime — 1 < A < 1. We 
give the free fields realizations of the boundary vacuum state and it's dual. Using 
these realizations, we give the integral representations for the correlation functions. 



^ ! 1 Introduction 

The one-dimentional massless spin | XXZ chain with a boundary is a system described 
by the Hamiltonian, 

- oo 

H * = -9 + a n+^n + A <+1<) + K, " K A < 1. (1.1) 

/ n=l 

Here the a®, o\ and a* stand for the Pauli matrices acting on the n-th site of the semi- 
infinite spin chain : 

••■<g>C 2 ®---®C 2 ®C 2 ®C 2 . (1.2) 
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In this paper we are interested in so-called massless regime : 

/ 7T \ 

— 1 < A = — cos I - — - J < 1, 4 : generic, 

where the spectrum of the Hamiltonian ( |1.1|) is gapless. In recent works |], [|, §, f| the 
various massless models "without boundary" were discussed, in the framework of the 
free field approach. In this paper we shall study the massless model "with a boundary" , 
in the framework of the free field approach. 

In the earlier work || the massive XXZ chain with a boundary was considered. The 
diagonalization of the Hamiltonian was obtained and the integral representations of the 
correlation functions were derived, in the framework of the free field approach. It's 
U q (sl n ) -generalization was achieved in ||. In the work 0, Baxter's Corner Transfer Ma- 
trix Method were extended to the boundary problem- the XYZ model with a boundary. 
The quantum Knizhnik-Zamolodchikov equation with boundary reflection, which governs 
the correlation functions, was derived. An infinite product formula of the one-point func- 
tion was derived by solving the difference equations. The Corner Trnsfer Matrix Method 
can be applied to only massive models. Fortunately, the massless XXZ spin is a limiting 
case of the massive model- the XYZ model fl|]. Therefore the correlation functions of 
the massless XXZ spin with a boundary were described by the following systems of the 
difference equations, which imply in particular the quantum Knizhnik-Zamolodchikov 
equation with reflections. 

G(Pi, ■ ■ ■ ,{3j+i,Pj, ■ ■ ■ , @2n)- ,e j+ i,£j,- 

R(P 3 -P^ie^G(P ir -- ,P v P j + ir -- ,P 2N )..., e ^ + i ,..., (1.3) 

and 

G((3i + 7tz, • • • , P2N-1, /32jv) ei! - ,e 2N = K((3i)Zl\G(-(3i + iri, f3 2 , ■ ■ ■ , A>iv)ei,- ,e 2N - 

(1.4) 

Here Gn{Pi, ■ • ■ , @2n) is a function with values in C® . The matrix R{(3) denots the 
/2-matrix defined in ( |2.1| ). The matrix K((3) denotes the boundary if-matrix defined in 
( p,14j ). In this paper we construct the integral reprersentations of the iV-point correlation 
functions, which satisfy the above systems of difference equations. 
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In this connection we should mention about the works B, |9|. They constructed 
the bosonizations of the boundary vacuum state of the type-II vertex opertaors, for 
the sine-Gordon model || and the 5'C/(2)-invariant massive thirring model ||, by using 
Lukyanov's bosonizations of the vertex operators with uv-cutoff [TO]. Therefore their 
constructions started with bosons with uv-cutoff, and the form factors were derived after 
removing the cutoff parameter at final stage. In this paper we prefer to work directly 
with operators with the cutoff parameter removed 0, and construct the bosonizations 
of the boundary state associated with the type-I vertex operators. 

Now a few words about the organization of the paper. In section 2 we formulate our 
problem. In section 3 we construct the bosonizations of the boundary vacuum state and 
it's dual state. In section 4 we derived the integral representations for the correlation 
functions. In Appendix A we summarize the bosonizations of the vertex operators Q. 
In Appendix B we summariz the Multi-Gamma functions. 



2 Formulation 



The purpose of this section is to formulate the problem. 
Let us set the i?-matrix as 



R{(3) = r(/3) 



HP) c(p) 

c(P) b{(3) 



(2.1) 



where we set the components as 



sh 



609) 







sh 



/3 + iri 
£ + 1 



Here we set 



sh 



sh 



£ + 1 



r(0) 



(2.2) 



(2.3) 



S 2 (i(3\27i, n(C + ljjjgHg + 7r|27r, 7r(e + 1)) 
~S 2 {-il3\2n, vr(e + l))S 2 {i(3 + n\2n, n(£ + 1)) ' 

where ^(/^iu^) is the double sine function defined in Appendix B. 

Let {v + , f_} denote the natural basis of V = C 2 . When viewed as an operator on V ® V, 
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the matrix elements of R(/3) are defined by 



R(P)v kl ®v k2 = E v h ® v h R {P) 
ji,h=± 

The i?-matrix satisfies the Yang-Baxter equation : 



k\k2 

hh ' 



(2.4) 



Rn(Pl ~ fo)Rl3(Pl ~ Ps)R2 3 (P2 - Ps) = R23(P2 - Ps)Ri 3 (Pi - (h)Ri2(Pi ~ (h) ■ (2.5) 
The normalization factor r (/3) is so chosen that the unitarity and crossing relations are 



R(-p) 



hh 



The commutation relation of the type-I vertex operator is given by 



fci,fc 2 =± 



The bulk scattering matrix is given by 

/ 1 

= s((3) 



\ 



V{0) c'{(3) 



where we set the components as 



sh 



in — (3 



sh 



in — P 



Here we set 



*(P) = 



i ) \ t 

S 2 (-iP\2n, n£)S 2 {n + i[3\2n, <) 



5 2 (i/3|27r,7rO-5 2 (7r-i/3|27r,7rO ' 
The commutation relation of the type-II vertex operator ^j(/3) is given by 

k\, k 2=± 

The commutation relation of thpe-I and type-II vertex operators is given by 



(2.6) 
(2.7) 

(2.8) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



The free field realization of the vertex operators on the Fock space of boson was given 
in I]. We summarized the bosonization of the vertex operators in Appendix A. 
Let us set the boundary i^-matrix by 

( 1 \ 



K{(3) = k((3) 



sh 







V 



sh 



where the normalization factor is given by 



(2.14) 



k{(3) = k (P)h(P), 



(2.15) 



where 

= S 2 (-2ip + tt(£ + l))S 2 {2ij3 + 37r|47r, 7r(^ + 1)) 

o[P) ~ S 2 (2if3 + 4vr|47r, vr(£ + l))S 2 (-2i/3 + 37r|4vr, vr(£ + 1)) ' 
S 2 (-i0 + + n\2n, Tt(£ + 1))S 2 (0 + in + 2n\27r, ?r(£ + 1)) 
X ^> ~ S 2 (if3 + i{i + 7r|27T, vr(£ + l))S 2 (-iP + ifi + 27r|27r, vr(£ + 1)) ' 

The matrix elements K{(5)^ are defined by 

j=± 

The i?-matrix and the .fT-matrix satisfy the Boundary Yang-Baxter equation. 



(2.16) 
(2.17) 

(2.18) 



K 2 (p 2 )R 21 (pi + i3 2 )K 1 {0 1 )R 12 {i3 1 - (3 2 ) = it! 21 (/3 a - faK^Rnfa + (3 2 )K 2 (f3 2 )(2.l9) 

The normalization factor k(f3) is so chosen that the boundary unitarity and the boundary 
crossing relations are 



K(f3)K(-f3) = id, 



k=± 



TV l 

~2 



Let us set the renormalized transfer matrix 



(2.20) 
(2.21) 



T B {p)=g-iY i *)(-P)KW i i*iW- 

j=± 



(2.22) 



Here we set the dual vertex operator as 



$*(£)= <M/3 + (j = ±). 
5 



(2.23) 



The constant factor is given by 



'Asin^) (, K+1) r( F l T )) 2 S tol, (,24) 



where g((3) is given in ( |3.18| ). 

The renormalized transfer matrix has the following relations. 



[T B (p 1 ),T B (p 2 )]=0, (ft,A6R), (2.25) 
T B (0) = zd, T B (p)T B ( r p)=id, (2.26) 
T B {-f3 + m)=T B {f3). (2.27) 



The Hamiltonian 7Yb ( |1 ■ 1| ) and the renormalized transfer matrix T B {j3) ( |2.22| ) are related 
through the formula, naively. 

d 



T B ) (0) ~ H B . (2-28) 



Inspite of constructing the eigenstate of the Hamiltonian ([1 .If) , we solve the following 
eigenstate problem of the transfer matrix T B (f3). 

T B (P)\B) = \B), (2.29) 
(B\T B (P) = (B\. (2.30) 

3 Boundary state 

In this section we invoke the bosonization method to find the explicit formulae for the 
boundary state \B), assuming uniqueness. The boundary state is determined by the 
following relation. 

T B (P)\B) = \B). (3.1) 
Actiong the vertex opertaors <f>j(—/3) from the left, we have the equivalent relation. 

K(pp j (p)\B) = *j(-fl\B), (J = ±)- ( 3 -2) 
Here we have used the duality relation. 

Qjifl&M = g x 6m, (j,k = ±), (3.3) 



where g is defined in Q2.2-J ). 



We make the ansatz that the boundary state has the following form. 

\B) = e F \vac), (3.4) 

where 



2J [b(t),b(-t)} 
When we set the coefficients as A(t) = — 1 and 



b{ - t)2M+ rmh\ b{ - t)dt (3 ' 5) 



m i « (?)*(«» -ft) i -(t)-(t)^K (36) 

B \ t ) = 7 Ttt \ 7 Ttt \ ' ( S - b ) 

1 sh (-(£ + !)*) * sh^+l)*) 



the boundary state |5) satisfy the characterizing relation 
Let us prove the relation (pT^). In what follows we use the abberiviations : uj\ = 2n, uj 2 = 
7r(£ + 1), and 

U + (,S) = exp (- jf 'M.e^dtY = exp (J™ ^e^'diV (3.7) 

°^ - or S-^) • °- {a) = exp (- r b -w e "" dt ) • <3 - 8) 

In what follows we omit non-essential constant factors. 
At first we explain the formulas of the form 

X{fc)Y{fa) = Cxyifc - fa) : X(fc)X{fa) :, (3.9) 

where X, Y — Uj, and Cxy{P) is a meromorphic function on C. These formulae follow 
from the commutation relation of the free bosons. When we compute the contraction of 
the basic operators, we often encounter an integral 

F(t)dt, (3.10) 



o 

which is divergent at t = 0. Here we adopt the following prescription for regularization 
: it should be understood as the countour integral, 

m^^dt, (3.11) 
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where the countour C is given by 




a 



Contour C 



The action of the basic operator on the boundary satate is given by 

U+(P)\B) = C(mnt.h(P)U-(-P)\B), (3.12) 

where 

h ,p. T 2 {-2i(3 + 47TI2W!, uo 2 )T 2 {-2i(3 + tt(£ + 1)|2^, u 2 ) 



T 2 (-2i(3 + 3tt|2u;i, cj 2 )r 2 (-2z/3 + 7r(^ + 1) + 7r|2u;i, w 2 ) 
x r 2 (-z/3 + tt|u;i, a; 2 )r 2 (— - t/j + ?r(£ + 1) + tt)^, u; 2 ) 
r 2 (-i/3 + i/i + 27r|o;i, uo 2 )T 2 (-i(3 - i/i + 7r(£ + l)|o»i, cj 2 ) 

The function satisfies 

= fc(/3) = W)*a(/3) = ^t^. (3.14) 

We have 

h{-p)$ + (J3)\B) = Const.h{-f3)h{f3)U_{f3)U_(-[3)\B). (3.15) 
Now we have proved the "+"-part of characterizing relation of the boundary state (fj.2|). 

K(p)X9+(p)\B)=*+{-P)\B)- ( 3 - 16 ) 



Next we shall prove the "—"-part of (|3.2|) . 



The commutation relation of the basic operator is given by 

U+WU-ifc) = g(Pi - /3 2 )C/_(/3 2 )C/ + (/3 1 ), (3.17) 

where we set 

a(B) _ r 7wArr ^Hff + 2 *\^ U2)r 2 (-t(3 + 7T(e + 1) K CJ 2 ) 

yW r 2 (-^/5 + 7r|a; 1 ,a; 2 )r 2 (- ? /? + 7r(e + l)+7r|a; 1 ,o; 2 )' l ' } 
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The action of the basic operator on the boundary state is given by 



U+(a)\B) = Const.I(a)U4-a)\B), (3.19) 

where we set 

1(a) = - *° : [ T\ (3-20) 



TC(£ + 1) r[ ili-ia. j 1_ 

1 l^TT) + 2(e+ 



2(^+1) 

From direct calculation, we have 



m ~ lsh (jTi) 

Const, xshf^ fda TT ^ + C f + — ^ 



■p / —IjA—lQL | -i 

x sh — 7T r x a 



-ifi—ia I -I 1 



x C/_(/5)[/_(-^)C/_(a)C7_(-a)|5). (3.21) 

Note that the operator part LL(/3)LL(— /3)C/_(a)C7_(— a) is invariant under the change 
of variables a <-> —a, /? <-> — /?. 
We get 

M/r 1 ^ f |±|) $-(/3)is) - M-zr^h fef ) $-(-/3)is> 
U + i/i-oo J 2 i ± v »r(e+i) 2(e+i)>/ 

TTrf ^ + ^ +l- * x aTTsh 
11 \ tt(p -u -h 9#-t 11 



id" V tt(£+i) 2(e+i); --ii""ve+i 2(e+i) 

x t/_09)C/_(-^)tL(a)CL(-a)|S). (3.22) 

The integrand of (RHS) of the above equation is anti-symmetric to a change of the vari- 
able a «-> —a. It means the left-hand side becomes zero after taking integral. Therefore 
we arrive at the following. 

k-p)& (||y) M/3)|5) = M/5)sh (frf ) ^ 

Now we have proved the "—"-part of the characterizing relation of the boundary vacuum 
state Q . 
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The dual boundary state (B\ is determined by the following relation. 

(B\T B (P) = (B\. (3.24) 
Acting the dual veterx operator $*(/3) from right, we have the equivalent relation. 

+ m)K{ 1) < = (B\^(P + ttz), (j = ±). (3.25) 
We make the ansatz that the boundary state has the following form. 

(B\ = (vac\e G , (3.26) 

where 

G=- [°° - C ^ b(t) 2 dt+ r rwfff m ( 3 -27) 
2 7 [6(f), &(-*)] 1 ' 7o P»W.6(-*)1 

When we set the coefficients as C(f) = — e~~ 2lTt and 

^ KD-'O^f-^ .^ ^^^T)-'^) (328 , 

* sh(-(^ + i)t) * sh g(e+i)t) 



the state (_B| satisfies the characterizing relation ( |3.25| ). It can be shown as the same 



manner as the case of the state \B). Here we omit details. 

4 Correlation functions 

In this section we calculate the vacuum expectation values of type-I vertex operators, and 
obtain them as integrals of meromorphic functions involving Multi-Gamma functions. 
We shall consider the 2iV-point functions defined by 

r fa a 

Lr ei ... e2N {P 1 , • ■ ■ ,P2N) — (B\B) ' 

From the commutation relation of the vertex operators fl2.8|) , the vacuum expectation 



values ([4 .1]) satisfy the i?-matrix symmetry ( |1.3|) . From the reflection relations fl3.2|) and 



( |3.25| ), the vacuum expectation values (|4.1p satisfy the reflection conditions 

Specializing the spectral parameters, they give multi-point correlation functions of 
the local spin operators of the massless XXZ spin with a boundary : 

G- ei ,...,- £jv , eiV ,..., ei (A + vrz, • ■ ■ ,0 N + iri,/3 N , ■ ■ ■ (4.2) 
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After specializing the spectral parameters as the above, our integral representations of 
the correlation functions ( |4.1|) can be compared with the formulae for the correlation 
functions of the boundary XYZ spin, which can be derived by mapping to the boundary 
SOS model |TTJ. Y. Hara |T^] considered the mapping method |]T3| of the boundary 



XYZ model, and derived the explicit formulae of the one-point functions, whose Ap- 
point generalization seems to be tedious but straightforward fl4]l . Comparsion with two 
formulae is our future problem. 

Now let us calculate the vacuum expectation value ft4.1| ), explicitly. 
Fixing indexes {ei, • • • , e 2 7v}, let us denote by A the index set 

A = {a\e a = -, 1 < a < 2N}. (4.3) 

In order to evaluate the expectation vbalue (3~I), we invoke the bosonization formulae of 



the vertex operators and the boundary state. By normal-ordering the product of vertex 
operators, we have the following formula. 



n 

l<bi<62<2Af 

TT r da nrf ^"^ i 1 W ^- a °) i 1 

XX / aa a XX i (f . u + 1 J M / 1 _i_ i \ + 



r 2 (»(Aa - 


AJ 


+ 27r|u;iu; 2 ) T 2 (i(A 2 


-AJ 


+ 7r(£ + 1) Co>iCU 2 ) 


T 2 (i(/3b 2 - 




) + tt|u;iu; 2 ) T 2 (i(A 2 


-Ax) 


+ 7r(e + 2)|wio; 2 ) 



X 



7r(e+i) 2(e+i); we+i) 2(e+r 

r 



n \ ^ ~ 



7T 



(e+i) e+i 



„! „.. I I I H aa 2 Qa l) I 



7T(e + 1) ^ + 1 

7r(e+i) ' 2(^+1)7 tt ~ we+i) ' 2(e+i; 



x J({(3 b } 2 b ^\{a a } aeA ). 



" ' , ^(A ~ ««) j _ _J__] i L a r f ijag ~ A) j _ 1 

. e A,i< 6 <2iv ^(^ + 1) ^ 2 (e + i)y a6W - 2iV W£ + i) 2(^ + 1) 



(4.4) 

Here we set 

(£| exp ( / X(t)b(-t)dt) exp ( / Y(t)b(t)dt J |£) 

A{A}^liKW) = ^ ^ U ° L , (4-5) 
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where 



2iV 



-iPbt „—ia a t 



x{t) ■ S=w-S^- (4 ' 6) 



Y ^ = -EiK^y + STTM = " x * (t) - (4 - 7) 

6=1 Sh ^^ aeA Sh (~t) 

Next we evaluate the quantity J({/3 b }\{a a }). Using the completeness relation of the 
coherent states (fjj, and performing the integral calculations, we have 

J({Ml^\{0Ca}aeA) 

1 sh(f)sh(7rt)sh(^) 



exp 



o l-A(t)C(t) tsh(f(e + l)t) 



x Qc(f)X(f) 2 + A(f)C(t)X(f)V(f) + l -A(t)Y(tf 

+ J°° 1 _ {(£>(*) + C{t)B(t))X{t) + (B(t) + A(t)D(t))Y(t)} dtj (4.8) 

Here A(t), B(t),C(t) and -D(t) are the coefficient functions in the boundary state and 
it's dual. It is evaluated as follows. 

J({Pb}\{<*a}) = M{fo})J a ({<Xa})Jfia({Pb}\{aa}). (4.9) 

Here we set 



j (fg\\ TT S 2 (ifi + i(3b + Tr\uj 1 uj2) tt / S 2 (2if3 b + 4ir\2uj 1 uj2)S2(-2i(3 b + 2it\2uj 1 u;2) 
6=i ^(^ + ifa + 27t\uj 1 oj 2 ) y S 2 (2t(3 b + 3tt|2^ 1 o; 2 )5' 2 (-2z/3 6 + 7r|2^ 2 ) 



x 



X 



X 



2N 

n 

6=1 



S 3 (2if3 b + 3711^1^1^2) S 3 (-2i/3 b + 7r|co>iCo>iCo>2) 
S 3 (2ij3 b + ATT\iu 1 uj 1 uj 2 )S 3 (-2i/3 b + 27t|cg>iu;iCo>2) 



n 

l<6i<6 2 <2Af 



S 3 (i(/3 bl 


+ (3 b2 ) + 3tt\uiUiuj 2 ))S 3 { 




+ (3 b2 ) + TT \U1U1U2)) 


S 3 (i(f3 bl 


+ (3 b2 ) + A^uxUiU-i^S^ 




+ P b2 ) + 2vr|u; 1 u;iu;2)) 



n Trior- fa « _l 1 u 5 3MAi-4 +^k w i w 2 1 
,_^o^M I S 3 {ie{(3 bl - (3 b2 ) + 47r|^i^i^ 2 ) J 

(4.10) 



l<6i<6 2 <2AT e=± 
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6=1 aGA 
2N 



X 



X 



nnn — 1 — (4.11) 

2N 

nnn 



S 2 (ie(a a - fa) + f |wiw 2 )ri(ie(oi! - /3 b ) + \ \uj 2 ) ' 



and 



mm) = n 



sin I Iff 



X 



X 



X 



X 



aeA sin ^ ?+1 + 2(?+1) 

nn 



eA e= ± ^/S 2 (2iea a + it\ujiuj 2 ) a/ S 2 (2iea> a + 7r(£ + 1) 
1 1 



n 

a€A 



Ti(-2m + 7r£|cu 2 ) y/S 2 (2ia a + tt(£ + 1) + vr|f uj 2 ) ^S 2 (-2ia a + 



11 8111 — + 



a x <a 2 e =± 



n 



, + 1 e + 1 

i 



X 



X 



o AA S 2 (*(a ai + a a2 ) + ^\fu; 2 )S 2 (-i(a ai + aj + tt£ + 2n\fu 2 ) 

a l > a 2 ^ ^ 

„n n 52( , e(aai _ aa2)+7r(e+1) | >2) 

°S S ri ( ie ( a «i ~ O + + l)|w2)ri(ie(a ai - a„ 2 ) + 7t£|cj 2 ) ' ^' l2 ^ 

Here we omit an irrelevant constant. 

In order to get the integral representations of the form factors of the local spin 
operators, we have to calculate the vacuum expectation value of both type-I and type-II 
vertex operators, and obtain them as integrals of meromorphic functions of Multi-Gamma 
functions. 

Caluculation of the vacuum expectation values is tedious but straightforward. We can 
perform it as the same manner as the correlation functions. 
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A Vertex Operators 

Here we summarize the bosonizations of the vertex operators PJ . 
Let us set free bosons b(t)(t G R) which satisfy 

sh ( — ] sh(-7rt)sh — - 
W), 6(01 = + 2 *(* + 0- (A.1) 



Let us set a(t) by 



K0*^ = a(t )sh^. (A.2) 



Let us consider the Fock space TC generated by the vacuum \vac) which satisfies 

b(t)\vac) =0 if t > 0. (A.3) 
The bosonization of the type-I vertex operators is given by 

d> + (/3) = U({3), (A.4) 



$_(/3) = J da: U(0)U(a) : 

Vtt^ + I) 2{£+l)J V TTf+l 2(^ + 1)7 ' 1 ; 
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where we have set 



U(a) =: exp (- £ , £,(«) =: exp (£ | e "*) : . (A.6) 

The bosonization of the type-II vertex operators is given by 

<M0) = (A.7) 



#_(/?) = J da: V(P)V(a) : 



where we have set 

= : 0- sS e '° Vt J : ' K(a) : exp {- L W" dt ) ' ' (A ' 9) 

Here the integration contours are chosen as follows. The contour Cj is (—00,00). The 
poles 

a - (3 = — + mr{£ + (n e N) (A.10) 

of r + 2(|+I)) are above (7/ and the poles 

a -/3 = nir(£ + l)i, (n e N) (A.ll) 

2 

of T ^— ^(g+i) + 2(|+T) ) are below Cj. The contour Cu is (—00, 00) except that the poles 

a-/? = +mr&, (n e N) (A.12) 

2 



of T — ^) are above Cjj and the poles 



7T2 

a-P= — -nn£i, (n e N) (A.13) 



of T (-*^=£1 - 4) are below C/j. 
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B Multi Gamma functions 



Here we summarize the multiple gamma and the multiple sine functions. 
Let us set the functions Fi(x\ui), r 2 (a:|a;i, 002) and T 3 (x\ui, 002,003) by 

K*r.(*M + 7B..(*M = j c ^-"^, (B.i) 

logr 2 (, kl ,. 2 ) - lB 22(xM = j c ^-r \l-e^l- e -^ Y (B ' 2) 

7 o ,_l \ f „-xt log(-t) 



7 /" (it 

logr 3 (x|wi, cj 2 , ^3) + ^B 33 (x\lo 1 , u 2 ,uj 3 ) = / — -e 



2nit (1 - e-^*)(l - e-^*)(l - e~^t) ' 

(B.3) 

where the functions Bjj(x) are the multiple Bernoulli polynomials defined by 



n;=i(e^-i) 

more explicitly 



E—B r>n (x\uj 1 ---uj r ), (B.4) 



n=0 



fln(*M = (B.5) 

LO L 

x 2 ( 1 1 \ 1 1 ( U0\ L0 2 



B 22 (x\u) = _- _ + _ + ! . (B.6) 

001002 \oo\ lo 2 J 2 b \002 oj\ 

Here 7 is Euler's constant, 7 = lim„^ 00 (l + | + | + '"" + ^ — logn). 
Here the contor of integral is given by 




Contour C 
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Let us set 



We have 



Si(x\u) = : — — . . — r, (B.7) 

o f I x r 2 (o;i + uj 2 -x\u^U2) fu Q x 

S 2 {X\UJ 1 ,UJ2) = ^T~\ \ ' B - 8 

5 3 (x|wi,W 2 ,W3) = TTJ 1 w^TI \ ( B - 9 ) 

r 3 (cji + u 2 + to 3 - x\uji,uj 2 , u 3 )T 3 (x\ui, u 2 , uj 3 ) 



T 1 (x\u) = e^-^s-Ii^M S!(x\u) = 2sin(7rx/cu), (B.10) 

V27T 



r 2 ~t~ ^iI^i^ ^2) 1 s 2 (yX ~\- uj\ \ uj\ , 1 i^x 



r 2 (x|wi,w 2 ) Ti(a;|w 2 )' S2{x\u) 1 ,u 2 ) Si(x\u 2 y Y^x^) 



x(B.ll) 



r 3 (a: + ^1^1,^2,^3) _ 1 ^(x + Ux\ui, u 2 , u 3 ) _ 1 , . 

r 3 (x|o;i, uj 2 ,uj 3 ) r 2 (x|a;2,a; 3 )' S 3 (x\u 1 ,u 2 ,u 3 ) S 2 (x\u 2 , u 3 ) ' 



\ogS 2 {x\^ 2 ) = J c -A^^^-]og{-t) — , (0 < Rex <co 1 + uj 2 ). (B.13) 



2-7T 

S 2 (x\u 1 u 2 ) = -^=x + 0(x 2 ), (x->0). (B.14) 

y/UJlU} 2 



7TX 7TX 

S 2 (x\(jJiuj 2 )S 2 (—x\ujiuj 2 ) = —Asm — sin — . (B.15) 

OJl uj 2 
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